Helium- and Lithium-like ionic sequences: Critical charges 
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Abstract 

In non-relativistic quantum mechanics we study the Coulomb systems of infinitely massive center 
of charge Z and two-three electrons: (Z,e,e) and (Z, e, e, e). It is shown that in both cases the 
total energy curve in Z is smooth, without any visible irregularities. Thus, for both systems the 
physical integer charges Z = 1,2, ... do not play a distinguished role as would be associated with 
charge quantization. By definition, a critical charge Z cr is a charge which separates a domain of 
the existence of bound states from a domain of unbound ones (continuum). For both systems the 
critical charges are found, Z cr ^ e = 0.91085 and Z cr ^ e = 2.009, respectively. Based on numerical 
analysis, the Puiseux expansion in fractional powers of (Z — Z cr ) is constructed for both systems. 
Our results indicate the existence of a square-root branch point singularity at Z cr with exponent 
3/2. A connection between the critical charge and the radius of convergence of 1/Z-expansion is 
briefly discussed. 
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Let us consider the Hamiltonian which describes the charged center Z and k electrons 

a=l a=l a<b 

The case Z = k corresponds to a neutral atom and Z = k + 1, 2, . . . describes positive ions. 
Sometimes, there exist negative ions Z < k, e.g. at Z = 1,/c = 2. Critical charge Z cr 
is a value of Z which separates the domains "existence/non-existence" of square- integrable 
solution(s) of the Schroedinger equation, T-C^> = E(Z)^>. In other words, it separates the 
domain, of charge Z, of the Hamiltonian where the bound state(s) exists from the domain 
where it does not. Constructively, Z cr is determined by a condition of vanishing the ioniza- 
tion energy. The goal of this short Note is to find the critical charge Z cr for two electron 
(Helium-like) and for three-electron (Lithium-like) sequence. It was claimed long ago that 
the question about critical charge is closely related to one of the first theoretical questions 
of newly born atomic physics about the radius of convergence in the 1/Z-expansion. 
Making change of variables in ((T]), r — > r/Z, we arrive at 

^ = -*X> - t- + iE- (2) 

2 f-f ^r a Z rob 

a=l o=l a<b 

where the new energy E(£ = j?) = ^p-- In general, for all studied cases the ground state 
energy E(£) is a smooth, slow-changing real function in £ G [0,£ cr ]. Since time of Hylleraas 
it was said that for two-electron case, k = 2, the radius of convergence for the expansion 
of E, 

oo 

E = Ys e "C (3) 

n=0 

coincides with = 1/Z cr . Furthermore, it was claimed that the nearest singularity to £ = 0, 
which defines the radius of convergence, is situated at real £-axis and even at £ = = £ cr = 
1/Z cr (see [l| and references therein where extensive discussion with extended bibliography 
together with historical account are presented). Based on the analysis of large number of 
computed coefficients of 1 /Z-expansion (~ 401) it was found that the critical charge for 
two-electron sequence is 

Z% e = 0.91103... , (4) 

and the singularity associated with it is an essential singularity of quite complicated nature 
[l|. Although this value of the critical charge is actually quite close to the first estimate of 



the critical charge given by Stillinger and Stillinger 

Z% = 0.8941... , (5) 



it is in contradiction with the nature of the singularity. It is claimed in |2j (see Eq. (2.17V 
that it is a branch point with exponent 3/2. Recently, the result (J3J) was challenged in [3] 
by using the accurate numerical analysis of a different set of coefficients in 1 /Z-expansion 
than one used in It is stated in [3] that 



Z% e = 0.90223... . (6) 

J] on the nature of the singularity as a branch point with exponent 
3/2. Our analysis unrelated to perturbation theory excludes the results flS])-® about a 
value of a critical charge as well as about a position of the singularity defining the radius 
of convergence of (j3]j. Our critical charge is close to z|r (see below). However, we confirm 
the observation by pj, {^J that the singularity at the critical charge is a branch point with 
exponent 3/2. We are not certain about a nature of singularity (ies) which define the radius 
of convergence of the 1/Z— expansion ([3]). 

Our approach is not related with an analysis of 1/Z-expansion being instead based on 
accurate calculation of the total energy as a function of the charge Z in vicinity of the critical 
point with further extrapolation to the critical point. Extrapolating function is assumed to 
be in a form of the Puiseux expansion for Z > Z cr , 

00 

E{Z) = J2 B n(Z-Z cr ) a " , (7) 

n=0 

with a condition that exponents a n should grow with the increase of n, a n < a n +i- Both 
cases of two and three-electron systems are considered. 

Helium-like sequence. In order to find the total energy of the l 1 ^ state (the ground 



state) as the function of Z we used as tria 



explicitly correlated functions (see e.g. {4]) 

N 

Atrial = 2^ Ai t exp (~ a * ri ~ ^ r2 ) + (1 2 )] exp (-0*7-12) , (8) 

i=0 

where {Ai} and {aij, a{\ are linear and non-linear parameters, respectively. It was shown 
that this basis provides at present the fastest convergence in ground state energy for Helium 



function a linear superposition of the exponential, 



{Z = 2) and H _ (Z = 1) among known bases. Namely using this basis the ground state 
energy for Z — 1, 2 was found with 24 s.d. {4]. 

Using the computer code kindly provided by Prof. K. Pachucki (see 
www.fuw.edu.pl/ krp/) designed following the work [4J we see that the fast conver- 
gence in energy is also obtained for Z other than 1 and 2. We calculated the ground state 
energies for different Z varying from 0.95 up to 2.00 with interval 0.05 with 12 significant 
digits (see Table I for some examples). It is worth noting that for Z < 0.95 the rate of 
convergence starts to deteriorate dramatically with a decrease of Z - it requires many more 
terms in OH]) to be included to reach desired accuracy unlike for Z > 0.95. We exclude the 
domain Z < 0.95 from calculation of energies based on ([H]). 

Calculations of the ground state energy for different values of £ (in 30 points) confirm that 
the energy i?2e(0 is a smooth, slow-changing function without any visible irregularities in 
domain [0, 1/Z cr i2e ]. It does not exhibit any distinguished role of points 1/M, M = 1,2,..., 
which would indicate the appearance of the charge quantization. 



TABLE I: Ground state energy E<i e for two-electron sequence for selected values of Z found using 

(fit) 



where all digits assumed to be correct and E)> ' from the fit (jiop . 



z 






1.3 


-1.029896662309 


-1.029896664 


1.25 


-0.933575272295 


-0.933575273 


1.15 


-0.756014315641 


-0.756014316 


0.95 


-0.462124684391 


-0.462124684 



We make fit with terminated expansion (J7J) whose length is subsequently increased from 
n = 3 up to n = 7. Taking a large number of points in the interval Z G [0.95, 1.5] and 
making fit it is obtained in a stable way with sufficiently high accuracy that a n takes either 
integer or half-integer values(!), 



«o = , «2 = 1 , otz = 3/2 , «4 = 2 , a 5 = 5/2 , a 6 = 3 , a-j = 7/2 



(9) 



It confirms, in particular, the rigorous result by Simon 5j about the absence of the square- 
root term, a = 1/2 and presence of the linear term a = 1 in this expansion. It contradicts 



a statement from 



6| that the term a = 3/2 is absent. In general, such an expansion is in 
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agreement with one proposed in 

as 

. Eventually, assuming the exponents fl9]) and taking 
data set for energy E 2e at nine points Z G [0.95, 1., 1.05, 1.1, 1.15, 1.2, 1.25, 1.3, 1.35] a careful 
interpolation leads to the expansion (see ([7])), 

Z 2 

E { l lt) {Z) = — ^ - 1.142552(Z - Z cr ) - 0.174110(Z - Z cr ) 3/2 - 0.7700097(Z - Z cr ) 2 

- 0.1399230(Z - Z cr ) 5/2 + 0.0224694(Z - Z cr f + 0.0087298(Z - Z cr ) 7/2 , (10) 
with the critical charge 

Z cr ,2e = 0.91085 . (11) 

It can be seen that the expression (flO]) reproduces 8-9 s.d. in all numbers for energies 
included in Table I. 

The explicit knowledge of the first terms of the Puiseux expansion ({TO]) allows us to 
check whether the singularity at the critical charge defines the radius of convergence of the 
expansion (J3]). In order to do it we construct the Puiseux expansion of the function E 2e (0 
near the critical £ cr = 1/Z cr , 

jZitor - O 3 + i (S 7 + ^ ~ C 2 r f ) &r ~ 7/2 + • • • , (12) 

where -Bo,i,4 are related with coefficients in the expansion (0), (fTUl) . It leads to the following 
form of the coefficient e n in ([3]) at large n, 

_ , , n r(|) 1 



r(n + i)r(|-n) grh 



B 3 \ 5 _ L ( B 5 2 B 3 \ 35 _ 2 

5 3 + ( B 5 - i cr - j s^tcr + [B 7 + ^ cr_ 2 _ - Lr Y ) ( 5 -2n)(7-2n)^ + ' 



(13) 



which is a type of 1/n— expansion for e n . Now we can make a comparison of these coefficients 
with ones calculated in [l|, see Table HT1 One can see that for n > 20, the coefficients are 
becoming sufficiently close. However, a deviation in coefficients is growing with an increase 
of n. For n = 300, where 1/n-corrections to the asymptotic behavior of e n at n tending to 
infinity presumably can be neglected, the deviation reaches ~6 times when the coefficients 
are more or less of the same order of magnitude 10~ 20 . If, for a moment, we exclude a 
possibility that the large-order coefficients in [l| might be calculated incorrectly (see e.g. 
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13] 



3j) we arrive at the immediate conclusion that the singularity related with the critical 
charge (fTTj) defined by the Puiseux expansion f|T2|) can not explain the large-n behavior of 
the e— coefficients. Hence, there must exist other singularity(ies) on the circle of convergence 
which define the asymptotic behavior of e n in at large n. 

TABLE II: Comparison the a n -coefficients in the expansion ([3]) calculated in [lj] (rounded to 3-5 
s.d.), obtained in [3J and found by using the formula (|13p . 



n 


m 


® 


(USD 


20 


-0.76862e-5 


-0.76862e-5 


-0.71492e-5 


100 


-0.398e-10 




-0.689e-10 


200 


-0.301e-15 


-0.222e-15 


-1.065e-15 


300 


-0.522e-20 




-3.396e-20 



Indeed, such singularities might naturally exist. Certainly, there exist the square- root 
branch points, which appear in Landau-Zener theory of level quasi-crossings (see e.g. 
due to the level crossing (for discussion see e.g. the case of quartic anharmonic oscillator 
8|, |9|). The most natural candidate for such a quasi-crossing might be a pair of two complex- 
conjugated square-root branch points due to the crossing of spin-singlet l 1 ^ and 2 1 S' states, 
which have the closest energies at real Z when both states coexist. Following this idea we 
assume that the behavior of the energy is 



£ 2e (£) = VT( + «) 2 + b 2 (At + A 2 [(£ + a) 2 + b 2 } + ...) 

near square-root branch points at a ± ib. We were unable to find a, b, A\,Ai with r = 
(a 2 + 6 2 ) 1 / 2 < — ^ — which would reproduce the behavior of e n coefficients at n > 100 found 
in [l|. It seems that the only possibility which is left is to assume the existence of three 
(or more) singularities on the circle of convergence: one of them is of the critical charge 
and the others are pair(s) due to level crossing(s). This possibility looks quite exotic and 
thus unlikely. Before going to explore it we think that the independent calculation of the 
large-order coefficients e n in is highly needed (see below Note Added). 

Lithium-like sequence. In order to find the ground state energy of (Z, 3e)-system as 



the 
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unction of Z we use the variational methods with Hylleraas basis set as a trial function 



11] . The ground-state wave function \1/ is expressed as a linear combination of i/j, the 



6 



antisymmetrized product of and the spin function x, 



-0 = A[<f>(f[,T$ 7 r%)x\ 

Atrial = K 1 ^V3 3 r"X|r 2 n |exp(-a 4 r 1 -^r 2 -Qr 3 )] , (14) 

X = a(l)/3(2)a(3)-/3(l)a(2)a(3) 

where {aj,6j,Cj} are positive real parameters and ri\^ are non-negative integers [ljj]. Using 
this basis set with 15 variational (non-linear) parameters plus analytic evaluation of some 



Hylleraas integrals as well as recursion relations, Puchalski et al [11] have reached the im- 
pressive accuracy 10~ 14 in the ground state energy for both Z = 3 (Li-atom) and Z = 4 
(the Be + - ion) which is the highest at present. The results of these authors also indicate 
that the accuracy of the previous calculatio n | lOl| is overestimated and is limited to 10 s.d. 
Making use the computer code employed in [ll|], and kindly provided by Prof. K. Pachucki 
(see /www.fuw.edu.pl/ krp/) and modified by adding MINUIT minimization routine from 
CERN-LIB we are able to see that the fast convergence in energy is also obtained for Z 
other than 3 and 4. 

Calculations of the ground state energy for different values of £ confirm that the function 
E-3e(0 is a smooth, slow-changing function without any visible irregularities in domain 
[0, 1/Z cr .3 e ] (see below). It does not indicate any distinguished role of points 1/M,M = 
3,4, . . ., whose would be associated with the appearance of charge quantization. 

We calculated the total energies for different Z varying from 2.02 up to 3.00 and reached 7 
s.d. with total basis length equals to 502 



161 ] . To make fit we assume the same exponents (Q 
in the Puiseux expansion. Data set which was used for interpolation contains seven points 
Z e [2.02, 2.03, 2.07, 2.08, 2.10, 2.12, 2.16] (see Table III). We arrive at the interpolation in a 
form of the expansion (EJ) 

E { 3 f e u) (Z) = -2.934278-3.390491(Z-Z cr . 3e )-0.114813(Z-Z cri 3 e ) 3/2 -1.102097(Z-Z cri3e ) 2 , 

(15) 

(cf. (llOp ). where the critical charge 



J cr,3e 



2.009 . (16) 



In order to check consistency we found that (fT5l) at the critical point Z = Z cr ^ e reproduces 
the ground state energy of 2e-system E 2e in 6 s.d. In general, the expression ( ITS"]) reproduces 



6-7 s.d. in all numbers for energies included in Table III. It is worth noting that the critical 
charge fflBI) is inside the intervals for critical charges proposed in 13j. It seems the method of 
analysis used in [13| are quite rough, leading to not very precise results - a direct calculation 
of the ionization energy oc (E 3e (Z) — E 2e {Z)) immediately exclude essential parts of those 
intervals for critical charge. 



TABLE III: Ground state energy E% e for three-electron sequence for selected values of Z found 
using (|14|) where all digits assumed to be correct and E^J 1 ^ from the fit f 1 1 5 [) . 



z 




E (fit) 


2.16 


-3.47810826 


-3.47810790 


2.10 


-3.25509127 


-3.25509091 


2.075 


-3.16479824 


-3.1647979 


2.02 


-2.97184 


-2.97184 



As a conclusion, we state that based on extrapolation of highly accurate results for the 



ground state energy for (Z, 2e) and (Z, 3e) to the critical charge, the branchpoint occurs 
with the exponent 3/2. This is in agreement with a statement made in Q-jsl for (Z, 2e). 



Furthermore, it is in agreement with recent results for one-two electron molecular systems 

n 

(2Z, e), (3Z, e), (AZ, e) and (2Z, 2e), (3Z, 2e) [14|]. All these results indicate a universal nature 
of singularity at critical charge for Coulomb systems. So far, present authors are unable to 
give a physics explanation of this phenomenon. For all studied systems we did not see any 
indication of charge quantization. 
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Note Added. Following a suggestion of one of the referees we made a preliminary nu- 
merical study of energy behavior vs Z of the 2 1 5' state of two-electron system in domain 
Z G [1,2]. These calculations were carried out using a computer code kindly provided by 
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Prof. K. Pachucki. We found the critical charge Z cr2 J = 1.02 (cf.flBJ)). Energy curves in £ 
did not display in a clear way a behavior indicating quasi-crossings, see Fig. 1. It might be 
considered as a signal that their square-root branch points are situated far away from the 
real £ axis. A localization of these branch points requires a separate study which might be 
done elsewhere. A general situation with analytic structure of E(£) remains unclear. 
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